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Abstract 

We provide the exact time-dependent cosmological solutions in the Randall- 

Sundrum (RS) setup with bulk matter, which may be smoothly connected to 

the static RS metric. In the static limit of the extra dimension, the solutions 

are reduced to the standard Friedmann equations. In view of our solutions, 

we also propose an explanation for how the extra dimension is stabilized in 

spite of a flat modulus potential at the classical level. 
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As a possible solution of the gauge hierarchy problem, Randall and Sundrum (RS) pro- 
posed an /Z2 orbifold model with non-factorizable geometry of space-time |1|], which has 
immediately attracted a great deal of attention. The model employs two branes, Brane 1 
(Bl) with a positive cosmological constant(or brane tension) Ai = QkilVP and Brane 2 (B2) 
with a negative cosmological constant A2 = 6k2M^, and introduces a negative bulk cosmo- 
logical constant A;, = —Gk'^M^. Bl is interpreted as the hidden brane and B2 is identified 
with the visible brane. Then the metric has an exponential warp factor which could be used 
to understand the huge gap between the Planck and eletroweak scales. Although the RS 
setup introduces cosmological constants k in the bulk and ki and ^2 on the branes, it still 
describes a static universe because of the fine-tuning between the bulk and brane cosmo- 
logical constants k = ki = ~k2, which is a consistency condition in the model. Hence, if 
the fine-tuning is not exact, the solution has the time dependence and the universe expands 
exponentially but its form is not suitable for the standard Big Bang universe after the 
inflation. To circumvent this cosmological problem, some approximation schemes regarding 
the brane matter in the static limit has been taken into account and(or) some conditions 
(such as the positive brane tension) for the brane and bulk cosmological constants are re- 
quired 1^. Thus Bl was thought of as the visible brane even though it is inconsistent 
with the original motivation for the gauge hierarchy solution ||l|. This problem, however, can 
be resolved with the Gauss-Bonnet interaction 0], bulk matter effects and extra dimension 
stabilization process 0, etc. 

In this paper, we will present exact cosmological solutions in the RS setup with bulk 
matter. Although there have been many cosmological solutions in the RS setup 0-^, the 
graceful exit problem from the inflation phase to the standard Big Bang cosmology has not 
been seriously considered yet. In addition, the role of the extra dimension in the presence of 
the bulk matter is not well understood. Our exact solutions converge to the RS metric if the 
space time is made to be static, and leads to the standard Friedmann equations if the fifth 
dimension is stabilized. In view of our exact solutions, we can find a clue for a stabilization 
mechanism of the fifth dimension and obtain a small compactified fifth dimension naturally. 
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Throughout this paper we consider a (4+1) dimensional universe with coordinate indexed 
by (0, 1, 2, 3, 5). The action describing the bulk matter as well as bulk gravity and brane 
matter is 

s^l d^x^g (f - A. + /:(^)) + E / ^'^f^^ - A.) , (1) 

;?=1,2 branes 

where we set the fundamental scale M = 1. and c!"^^ represent matter contributions 

in the bulk and on the branes. For compatibility with the cosmological principle that our 
three dimensional space is homogeneous and isotropic, we assume that the metric of the 
universe has the following form, 



(2) 



where r denotes time and y denotes the fifth component. Prom the metric ansatz, the 
Einstein tesor Gmn is derived through the standard calculation. 
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(8) 
(9) 



where dot and prime denote the derivatives with respect to r and y, respectively, and the 
i runs through 1, 2, and 3. Here diagonal Einstein tensors are split into two parts, G^^^ aa 
and G*^^-* aa, depending on the nontrivial y and r derivatives, respectively. Thus the original 
Einstein tensor is, of course, expressed as the sum, Gaa = G^^^ aa + G^"^^ aa, where the A is 
(0, i, 5). The source part of the Einstein equation is composed of the cosmological constant 
and the energy-momentum tensor of matter. In this paper, we will regard the matter as 
perfect fluid. For the future convenience, we divide also the source tensor into two parts. 



T(i) ^B = -(l-r,)- diag [A,, A^, A;,, A^, A;, 



- E ^{y-y^y ""diagiA^A^A^Aj-.O] 

j=l,2 branes 

+ E ^(y - 2/j)e"^diag[-pj,Pj,Pj,Pj-,0] 

J =1,2 branes 

T(2) ^ 5 = diag [-(p + r/Afc), P - r]A,, P - r]A,, P - r]A,, P, - r]A,] 



(10) 



(11) 



where the pj and pj are nontrivial components of the energy-momentum tensor of the matter 
hving only on the j-th brane, and ?7 is a number representing how A^ is spht into T*^^^ b 
and T(2)^5. The total source tensor is described as b = T^^^ b + T"*-^-* ^ b- Here we set 
7o5 = because it is believed that there is no flow of matter along the fifth direction. The 
continuity equation of the energy-momentum tensor b -A = ^ must be satisfied, whose 
5 = and B = 5 components are 



p + 3i (p + P) + 5 (p + P5) = 
P^ + 3A' {P,-P) + N' (P + P5) =0 



(12) 
(13) 



The B = i component is identically zero. 
Now let us take some ansatze, 



G(^) AA = AA (or AA = aa] 
A'iT,y)=N'iT,y) . 



(14) 
(15) 



The above ansatze have been chosen to fulfill our purpose of restoring the Randall- Sundrum 
metric in the static limit. The ansatz Eq. (p!4D and Cos read 
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A" + 2A'^ - A'B' 

= -(1 - r])A, - e-^ [6iy) (A, + p,) + 5{y - 1/2) (A2 + P2)] 

2A" + 3A'2 + A^" + + 2A'N' - 2A'B' - P'A^'" 

= -(1 - r/)A, - e-^ [6{y) {A, - p,) + 6{y - 1/2) (A^ - P2)] 



A'^ + A'N' 



il-v)A, 



A'B + An' -A'-Aa' = o 



(16) 

(17) 

(18) 
(19) 



Under the ansatz Eq. (|T5D, Eq. (pTSD becomes 

A'2^_(l_^)^xe2^ = A;V^ . (20) 

Hence, de Sitter(A5 > 0,?7 > 1), anti-de Sitter (A^ < 0,?7 < 1), and flat Minkowski space 
(Af, = 0) are possible in the bulk. The solution consistent with the jZ^ orbifold symmetry 
is 

A(r, \y\) ^ kFir, \y\) + J(r) and F(r, \y\y = -e^^^^^^hgniy) , (21) 

where the sgn{y) is defined as sgn{y) = \y\' = 2[9{y) — 9{y — 1/2)] — 1. Then, because of 
the ansatz Eq. (pTS]), the exponential factor of the goo component in our metric tensor is 
written as 

N{t, \y\) = kF{T, \y\) + K{t) kF{T, \y\), (22) 

where K{t) is removed by the redefinition of time r in the second part of the above equation. 
Therefore, we ignore K{t) below. 

The above result Eq. (|2l|) leads to some useful relations, 

A" = -ke^B'sgn{y) - 2 [6{y) - 5{y - 1/2)] ke^ 

= A'B' - 2 [5iy) - 5{y - 1/2)] ke^ (23) 
i' = -ke^sgn{y)B = A'B . (24) 

Eq. implies that the N{t, \y\) should be stabilized if the B{t, \y\) can be stabilized 
somehow, since A' = N'. With Eqs. (|l^) and (|2^), we can show directly that our ansatz 
is consistent with Eq. ([I9|) . Because of Eqs. (|l^) and (|23|) , Eqs. (^) and ( |17|) just require 
matching the boundary conditions, 

/c = ^ (Ai + pi) = -i (A2 + P2) and pj = -pj . (25) 

Hence, considering the fiuid continuity equation on the branes, pj + SA{pj +Pj) = 0, we can 
arrive at a result pj = constant and so the pj is a constant also, which are expected results 
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from our assumption T05 = 0. The non-dynamical properties of the pj and pj obstruct 
our intention to explain the dynamics of the space-time exactly with brane matter, which 
is a reason for us to consider bulk matter. [The solution with the brane matter condition 
Pj = —pj is equivalent to that we do not have a solution with the radiation or matter 
dominated phases by brane matter. Anyway, the radiation or matter dominated universe 
must results with a vanishing effective cosmological constant, in which case we have not 
obtained a most general solution.] 

Now that we have fulfilled the ansatz Eq. ([l^) already, the remaining equations, G*^^-* aa = 
T(2) ^^are 



p + riAh = 3e-^^ A' + AB (26) 



P-7]Ab = -e"^^ 2A + 3A^ + B + B'^ + 2AB - 2AN - BN (27) 



P5 - r]Ab = ~3e~'^ [A + 2A' - AN\ , (28) 

which describe the relation between matter and geometry dynamics. They are nothing but 
the extended Friedmann equations. With Eqs. (|T^) and (|2^), we can check that the above 
equations Eqs. (^), ( P7| ) and (^E^) satisfy both fluid continuity equations, Eqs. (|T^) and (|T3p 
identically, that is, the constraints, Eqs. (|7^ and ( |73| j are just redundant equations, which 
are interesting results.^ Therefore, the remaining required conditions for the solution are 
only Eqs. (|2T|) and (^21). The relations among the p, P and P5 may be, of course, governed 
by particle physics. 

Toward a simple solution, let us consider the case that the size of the extra dimension is 
stabilized, i.e. B = 0, which leads also to B' = generically by redefinition of y. Because 



of Eqs. (|T5D, ( p2D and (p^), then, N is also generically set to zero. After all we have 

B = B' = N = . (29) 



Compare with the standard cosmology, where two Freedmann equations and one continuity 
equation lead to one dependent equation. 
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Then, from Eqs. (pTD and (p2D , the function -F(r, |?/|) is determined to -F(t, |?/|) = —ke^\y\ 
and so 



Ar(r, |y|) = -A;e^|y| = -A;6o|l/| 
A{t, \y\) = -kbo\y\ + J{t) = -kbo\y\ + I H{t)dt , 



(30) 
(31) 

where the interval scale Bq is a small constant and H{t) is a time dependent arbitrary 
function but may be determined by the equation of state. Thus the metric is read off as 



-2kbo\y\ 



(32) 



Note that the metric is the same as that of the Randall-Sundrum except for the factor 
^2fH{t)dt_ rjnj^g^ ^ become 



-_,2kbo\y\ 



P{T,\y\)-r]kk = -e 

P,{T,\y\)-r^K, = -?,e^^'''\y 
1 



2H[t) + ?,H\t 
H{t) + 2/f2(r 
p(r, - 3P(r, lyl) 



(33) 
(34) 

(35) 



2 " ' 



where /i runs through 0, 1, 2, 3. [-B(r, |y|) is associated with the vacuum expectation value 
of a massless four- dimensional scalar field.] The above equations, Eqs. (0), (p4D and (pSj), 
show that due to the exponential factor matter in the hulk is accumulated mainly near the 
B2 brane (negative tension brane), which is a similar result to those of the Refs. which 
are the exact solutions of bulk scalar, gauge boson and fermion fields to their field equations 
under the fixed RS static background geometry. 

Of course, any H{t) with H{t) —>■ and H{t) ^0 as r ^ oo can lead to an exit from 
an inflationary phase to a static Randall-Sundrum (k 0) or Minkowski (k = 0) universe. 
In this paper, however, we will not specify a model because we are more interested in the 
real expanding universe. 
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In Eqs.(|33|), (|3^) and (|35|), we should remember that p + rjK},, P — ?7A6 and P5 — ?7Ab are 
non-trivial components of the 5 dimensional energy-momentum tensor. To derive effective 4 
dimensional energy- momentum tensor b, it is necessary to consider the definition of the 
original 5 dimensional energy-momentum tensor, 



<5^matter = / d'x^^ Sig"" a)T'^'^^ B = j d'x j dy 6ov^ 5(5^ a)T(2)^ 



B 



(36) 



where = det[g^^] (/i, u = 0,1, 2, 3). As the 4 dimensional metric at a 4 dimensional slice 
in the bulk is g^u = e^'^'^^^^^g^p, which was introduced by Randall and Sundrum to solve the 
gauge hierarchy problem the effective 4 dimensional energy- momentum tensor is given 
as 



/ dy e-^fc^-ol^^lT^^)^, 
6n / dy e-2'=^«lfl 



• diag[3i/'(r), 2i/(r) + ?>H\t), 2H{t) + ?>H\t), 2H{t) + 3H\r)] 
diag[-(p(r) + r]A), P{r) - r]A, P{r) - r]~A, P{r) - r/A] . 



(37) 
(38) 



According to RS, 60 / dye ^^^oIj/I jg nothing but the induced 4 dimensional Planck scale Mpi 



II]. Thus, from Eqs.(|3D, (H, 

p(^> \y\) + v^b 



3 

g-2fc6ol?/l 



and (|38|) , we can get the Friedmann equations, 
f 



3M2^ 



p(r) + r^A 



and 



p{T,\y\) + 3P{T,\y\)-2r]At 



1 



p(r) + 3P(r) - 2r/A 



(39) 



, (40) 



a(r) 6 

where a(r) is a scale factor of our three dimensional space, a(r) = e-f ^(*)'^*. Therefore, we 
could have saved the whole standard cosmological scenario by only requiring stabilization of 
B in our framework (i.e. B = 0). Thus, we do not have to require additional conditions 
except for 5 = such as the bulk to have non-zero cosmological constant or the visible 
brane to have a positive tension, which were essential in other models 

As an example, let us consider the vacuum dominated era. The equation of state is 
P = —p and then P5 is given by P5 = —2p — r]Ah. Then Eqs. (|33|) and (0) lead to 
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H = constant = Hq , (41) 

which gives the inflationary universe. Hq can be considered as a parameter representing the 
degree of fine-tuning as given in Eq. (^). If Hq vanishes, the fine-tuning is successful and 
the universe is static. On the other hand, a non-zero Hq does not satisfy the fine-tuning 
condition and gives rise to an infiationary universe. Note that our solution is for B = and 
any modulus potential is not generated at the classical level since 60 is arbitrary. 

We proceed to discuss the possibility of stabilizing the size of the fifth dimension in our 
framework. Since we obtained already the solutions for the B = case, we conclude that if 
B goes to zero asymptotically, there exist solutions converging asymptotically to our B = 
solutions. 

We suppose that the stabilization era occurs before or during the conventional infiation 
era. Thus, we suppose that the system is governed by the same equations of state as those 
of the infiation era discussed above. Then, from Eqs. (0), (0) and (|28|), we obtain 



2A + B + B^-AB~ 2AN -BN = (42) 
A - AN-2Ab = , (43) 



which are the equations of state in the infiationary era. The above equations and Eqs. (|2lD 
and (p2|) convince us that any potential for the modulus field is not generated also since the 
B is not fixed yet. 

Eq. (|43| ) is easily solved, 

A = e^i\y\)+N^2 ^ ^si\y\)+kF^2 f j^p ^ j\ (44) 



where s{\y\) is an integration constant and b is defined as 6(r, \y\) = e'^^^''^'-'. Note that A is 
an increasing function of time since A > 0. Removing A from Eqs. ( ^2]) and (|5), we obtain 

B + B^ + sAB -BN = 0, (45) 

which can be solved to give 
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where t{\y\) is an integration constant. Note that b is an increasing function of time but b 
could be zero asymptotically. The extra dimension scale b can be stabilized if the combina- 
tion — + 3A (= 2kF + 3 J) is an increasing function of time without limit. Especially, if 
N A and A increases as r ^ oo without a limit at least with a power law, which inflates 
the three dimensional space, then b would decrease exponentially and so b could be stabilized. 
Below we show that it is possible. 

From Eqs. and (0), we obtain 

^ ^ gS(|y|)+t(|s/|)+3A ^ ^^Y^ 

b 

which is integrable. The solution is 

6(r, \y\f = u^{\y\) - e--(ls/l)-*(l?^l)-3^ , (48) 
where is a \y\ dependent arbitrary function. The remaining constraints to satisfy are 



only Eqs. and (p^), i.e. A' = N' = —kb{T, \y\)sgn(y), from which Eq. (^8]) can be 
written as 

3 X u\\y\)-b^ + "'^'^'^ + ^'^'^'^ = = ■ '^""^y^ • ^'^^^ 

We intend to make our solution become the inflationary solution asymptotically that was 
obtained before. Since —kbo\y\, A Hq, and 6 — > 6o as r — > cxd, let us take the 

integration constants in Eqs. and (EHI) as 



si\y\) = kbo\y\+\n Ho-2\n bo (50) 
u{\y\) = bo , (51) 

where bo and Hq were defined above. So far the solutions were exact. Note that any value 
for Hq is possible, which does not influence Eq. (|49|). Although Eq. P9| ) is difficult to solve 
exactly, we can argue that |A^| = k\F\ ^ {^(t) and Hqt} is sufficient to draw a meaningful 
conclusion. 
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J(r) and Hq can be chosen always such that |A^| = k\F\ ^ {^(t) and Hqt}. It is 
compatible with the phenomenological requirement of a large Hq so that our solutions are 
derived during the inflationary epoch. (Actually during the inflationary era, Hq^ is (a very 
large value) ~^ ~ 10~^^ sec.) In this early inflationary era, and b are regarded as being 
much smaller than Hqt, since the scale of the universe was small right after the Big-Bang. 
Then — + 3A increases with time and the extra dimension scale b is stabilized rapidly to 
bo. Once b is stabilized, the conditions Eq. (|^) become valid and A^ and A are forced to 
—kbo\y\ and Hq, respectively. 

To see our argument explicitly, let us take somewhat large J(r) assumption and solve 
Eqs. (0), ( ^6]) and (|^) approximately under the condition k\F\ <^ J. Our aim is to show 
A Hq, 6 — > 0, and b b^. Here, we set kF and b as 0(1) initially, and then J ^ 1. Then 
Eq. (IID gives 

bl - b%r, \y\) = ^e-''"\y\-^'^e'^' « 1 . (52) 
With Eqs. (^) and (|52| ) we are led to the following results, 

6'<1 or 6^6(r), (53) 

and we obtain 

F{T,\y\)^-bir)\y\ . (54) 

Here we must set t{\y\) ~ in view of Eq. (^61). Then, Eqs. ( ^4]) and (|46|) become 

J ^ Hoj^ and b ^ e"^-^ < 1, (55) 

which shows that during the period of the three space inflation it is quite difficult for b to 
be dynamical. It is interpreted as the stabilization of the extra dimension in spite of the flat 
potential for b. From Eq. (55) we can derive an expression for 6(t), 



4'" 



(60 - >>)- 

+ b„b + V 
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bp + 2b 
boVS 



^r, (56) 



or 

6(r) ^ 60 - e'=^^«" . (57) 

Here we can see that as r ^ 00, 6(r) grows to Bq and J tends to Hq asymptotically. In other 
words, to obtain an inflationary universe b should be stabilized to bo exponentially. Note that 
Eq. (|56D becomes an exact result provided the warp factor vanishes, which corresponds the 
cases of Afo = or = 1. 

If b is 0(1) but small right after the Big Bang, bo should be 0(1) but small. As the extra 
dimension gets stabilized {B = 0) soon after the beginning of the inflationary era, while 
the three dimensional space inflates (to e^°'^), b remains small (~ 1/M) and the universe is 
reduced effectively to 4-dimension. 

In a similar method, we can show that b is made asymptotically to zero also in the 
radiation {P = p/3, rj = 0) and matter dominated era {P = 0, rj = 0). However, as the 
initial condition for b is zero (through the above solution in the inflationary epoch), b should 
have been stabilized already and so Eq. (^) should have been valid since the beginning of 
the radiation dominated era. Therefore, the Friedmann equations Eqs. (|39D and (40) hold 
good in the radiation and matter dominated eras. 

In conclusion, we have provided exact cosmological solutions in the RS setup with bulk 
matter. In the static limit of all components of the metric, the solutions become the RS 
metric and in static limit of the extra dimension, they are reduced to the standard Friedmann 
equations, which implies that bulk matter is accumulated mainly near the negative tension 
brane (visible brane B2). In this case the modulus potential is not generated effectively at 
the classical level. With our solution, however, we have shown that the extra dimension 
could be stabilized (the B = solution) even if the modulus potensial is flat (60 is arbitrary) 
and it should be small since the three dimensional space inflates during the inflationary era. 
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